"""
Created on 2018/4/10 9:55 星期二
@author: Matt  zhuhan1401@126.com
Description: 与数学性质相关的工具箱
"""

import numpy as np

def isPrimes(x):
    a=np.arange(1,x)
    nMax=int(np.sqrt(len(a)))
    isPrime=np.ones(len(a),dtype=bool)
    isPrime[0]=False

    for i in range(2,nMax):
        if i in a[isPrime]:
            isPrime[(i**2-1)::i]=False
    print(a[isPrime])

isPrimes(101)